Dynamics of Low-Density Ultracold Rydberg Gases 



OO 
O 

o 

(N 



43 ■ 

i ' 

> ■ 

C ■ 
c3 ■ 

^ : 
cr 



> 



o 

OO 
O 



J. 0. Day, E. Brekke, and T. G. WalkeiB 

Department of Physics, University of Wisconsin-Madison, Madison, WI 53706 

(Dated: March 7, 2008) 

Population dynamics in weakly-excited clouds of ultracold 87 Rb Rydberg atoms were studied by 
means of trap loss, fluorescence detection, and state dependent stimulated emission. Rydberg atoms 
were excited to various nl Rydberg states via continuous two-photon excitation from a magneto- 
optical trap. A stimulated emission probe laser was then used to bring the Rydberg atoms down to 
the 6P3/2 state, allowing state-dependent detection of the Rydberg atoms. Measurements of trap 
loss and fluorescent emission reveal information about the evolution of the Rydberg populations. In 
particular, population in the initial Rydberg state quickly transfers to other Rydberg states by a 
non-collisional mechanism, likely superradiant emission. The trap-loss measurements are consistent 
with black-body ionization as the dominant loss mechanism. 

PACS numbers: 32.80.Ee,42.50.Nn,34.50.Cx 



In the past 10 years it has become possible to study 
the interesting properties of dense clouds of ultracold Ry- 
dberg atoms, usually produced by laser excitation from 
magneto-optical traps (MOTs) fil m. A variety of inter- 
esting collisional phenomena (STEME [E 0, B B 03 H3 
have been observed, including surprisingly rapid sponta- 
neous conversion of dense ultracold Rydbergclouds into 
plasmas 0,112, El- It is clear that resonant "M 3, H, Bl 
and near resonant energy transfer collisions [8|, |9|, LLP] play 
important roles in the population dynamics of clouds of 
ultracold Rydberg atoms. 

In addition, intriguing proposals to use the strong, 
long-range interactions between ultracold Rydberg atoms 
to perform conditional quantum manipulations [14| are 
driving a number of groups to study coherent manipu- 
lations with small numbers of atoms and small samples. 
In addition, the concept of Rydberg blockade [l5[ brings 
up the possibility of using spatially confined samples of 
atoms (lfjl [13] for mesoscopic quantum manipulations. 
Using such collective Rydberg excitations might allow 
novel quantum applications such as fast quantum gates 
15L [l8| , single atom and directed single photon sources 
19(, and fast quantum state detection [20(. 

The ultracold Rydberg cloud experiments, with many 
Rydberg atoms excited at once, and the quantum ma- 
nipulation experiments, with only one Rydberg atom ex- 
cited at a time, operate in very different regimes. Still, 
it is interesting to note that in both cases their success 
is reliant on the strong interatomic forces between Ry- 
dberg atoms. For quantum manipulations, one wants 
the Rydberg-Rydberg interactions to be so strong that 
multi-atom excitations are suppressed plj . A number 
of experiments [12, LH, LH, [H| report suppression signa- 
tures suggesting that this regime should be achievable. 
Most recently, collisional dephasing of coherent Rydberg 
Rabi flopping was observed for only two atoms confined 
in a dipole trap 26] , and dephasing of coherent Rydberg 
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excitation in a BEC was also observed [2 71]. 

The experiment reported in this paper studies a tran- 
sition regime between the two types of experiments. 
We weakly excite ultracold Rydberg atoms in an un- 
blockaded regime and study their evolution using three 
tools: trap loss, spontaneous emission, and stimulated 
emission. Using these diagnostics, we find that under our 
weak excitation conditions the probability is very small 
that the Rydberg atoms experience an inelastic collision 
sufficiently strong to leave the trap. They predominantly 
return to the atomic ground state by spontaneous emis- 
sion. Using the stimulated emission diagnostic, we also 
observe rapid population transfer out of the initial Ryd- 
berg state-the initial Rydberg state is depopulated on a 
time scale substantially shorter than expected for trans- 
fer from black-body radiation or single-atom spontaneous 
emission. The transfer rate slowly decreases with increas- 
ing principal quantum number, strongly suggesting that 
the transfer process is not collisional in nature. We build 
upon the recent observation of superradiance by Wang et 
al. [28| and argue that superradiance is the likely mecha- 
nism for this transfer process. Finally, we show that trap 
loss rates are consistent with expectations from black- 
body ionization. We conclude with implications of these 
results for other experiments. 



I. TRAP LOSS STUDIES 

In this section, we describe observations of loss of 
atoms from the MOT due to Rydberg excitation. Other 
groups [H, [2i| have reported observations of Rydberg 
trap loss under various conditions, without attempting 
to explain the mechanisms for trap loss, especially un- 
der the low excitation conditions reported here. We find 
that trap loss rates are much less than excitation rates, 
showing that, once excited, the Rydberg atoms primarily 
decay back to the ground state by spontaneous emission 
processes. 

We excite magneto-optically trapped 87 Rb nlj Ryd- 
berg atoms using continuous-wave 5Si/a — * 5P3/2 — ►nlj 
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FIG. 1: Diagram of the Rydberg excitation scheme. 



two-photon excitation. The 780 nm (20 mW) and 480 nm 
(10 mW) lasers used are tuned Si — 500 MHz above the 
intermediate 5P3/2 F=3 state, thus avoiding direct exci- 
tation of the 5P 3 / 2 state (Fig. [I}. Both excitation beams 
have spatial dimensions on the order of the ~ 1 mm size of 
the MOT. Typical two-photon excitation rates vary from 
10 to 100/s depending on n level. These rates are de- 
duced from calculated atomic matrix elements, measured 
laser intensities, and the observed transition line-widths, 
and are averaged over the atomic and laser spatial distri- 
butions. The effective excitation rate is |e2 | 2 /A, where 
(-2 = e r ei,/4:Si is the two-photon Rabi frequency, e r and 
£f, are the single photon Rabi frequencies for the 780 nm 
and 480 nm lasers, and A is the observed transition line- 
width which is typically 8 to 10 MHz for these scans as 
described below. This results in ~10 4 Rydberg atoms at 
a relatively low density of 10 7 /cm 3 . This density is an 
order of magnitude or more smaller than the densities at 
which ultracold plasmas are formed 

BO- By chan gmg 
the frequency of the 480 nm laser, we excite to Rydberg 
states at 28D, 43D, 58D, and 30S; in most cases the D- 
state excitations are to the J=5/2 state. These particular 
states are chosen to provide a range of different sponta- 
neous decay rates (which scale roughly as ~ n~ 3 ) and a 
wider range of Rydberg-Rydberg van der Waals interac- 
tions Cgi?" 6 , with C 6 (28,43,58) w (0.08,540,330) GHz 
/Ltm 6 [2l|. 

A very simple and effective way to observe excitation 
of Rydberg atoms in a MOT is via trap loss, with a re- 
sulting decrease in the observed fluorescence from the 
MOT. Processes such as black-body ionization, photo- 
ionization by the applied lasers, and inelastic collisions 
between Rydberg atoms and ground state atoms either 
ionize the atoms or give the atoms sufficient (> 10 K) 
kinetic energy to leave the trap even if they radiatively 
return to the ground state. Radiative processes such as 
photon emission or absorption, elastic or inelastic col- 
lisions with cold electrons, and near-resonant Rydberg- 
Rydberg collisions do not transfer enough kinetic energy 
to the atoms to cause trap loss. Trap loss is very sensi- 
tive; it is easy to observe a 0.1/s change in the trap loss 
rate from the MOT. With modest Rydberg excitation 
rates of 10/s, a 1% inelastic or photo-ionization channel 
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FIG. 2: Excitation beams are scanned over the Rydberg res- 
onance at the 28D state. The peak excitation rate is 110/s. 
The figure shows the loss rate from the trap induced by the 
excitation beams as well as the photon count rate from the 
spontaneous emission cascade. 

can easily be detected. 

Since the excitation lasers are sufficiently off-resonance 
to cause negligible perturbations to the MOT trapping 
and cooling processes, the loss due to production of Ryd- 
berg atoms affects the MOT only through a change in the 
ejection or loss rate from the trap. If the MOT lasers pro- 
duce a loading rate L and a loss rate Tq, the steady-state 
number of trapped atoms is N g o = L/Tq. Excitation of 
Rydberg atoms using lasers tuned to frequency v adds a 
new loss rate Y\ which then changes the number of atoms 
to N g (v). The loss rate can then be found from 

This is operationally simpler than taking separate MOT 
loading transients at each frequency and determining the 
loss rate from the time constant. We have checked that 
the two procedures give the same results. 

Figure [2] shows the results of sweeping the two-photon 
excitation frequency through a typical Rydberg reso- 
nance. Loss rates induced by the excitation lasers are 
around 0.2/s, a rate much smaller than the 110/s Ry- 
dberg excitation rate. We may infer from this that the 
probability of inelastic collisions occurring with sufficient 
energy transfer for the atoms to leave the trap is less than 
1/500. For the other states studied, we find 1/50 for 43D 
and 58D, and 1/1000 for the 30S state. 

Noting that the trap loss probabilities are significantly 
higher for 43D and 58D as compared to 28D, a logical hy- 
pothesis would be that inelastic collisions are responsible 
for the trap loss. To this end, we compared the trap loss 
rates from 41D and 43D, and found little difference. Since 
the 43D van der Waals interactions are anomalously large 
due to a near resonance in the 41P+45F and 43D+43D 
potentials @, HH, [3(3] , then very different results should 
be obtained when comparing 41D and 43D if collisions 
were responsible. Thus we again find no evidence that 
resonant energy transfer collisions cause significant trap 
loss. 
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The following estimate supports these conclusions. 
The inelastic collision rate should be a capture rate mul- 
tiplied by a probability of energy transfer. We estimate 
the capture rate as rjva, where a — ttRq, CqR^ 6 = kT, 
and r\v is the Rydberg atom flux. We then deduce a 
capture rate of 

r/va = r/vTT {C e /kT) 1/3 ~ 200/s (2) 

for an n = 43 density of rj = 10 7 /cm 3 . This rate, al- 
ready small, is further reduced by the energy transfer 
probability. For states like 43D, where the van der Waals 
interactions are repulsive at long range, there are no ther- 
mally accessible curve crossings and we expect the rate to 
be suppressed by at least a Boltzmann factor e - AE / kT ; 
which even for 43D with small AE is a factor of ~ 10~ 2 . 

The line- widths observed are significantly broader than 
the laser line- widths, and exceed what we expect from 
drifts of our reference cavity used to stabilize the 480 nm 
laser frequency. This is consistent with observations by 
others [11 [H|. 

Since the atoms are only rarely leaving the trap due 
to Rydberg-Rydberg collision or ionization, they must 
primarily return to the ground state by emission of one 
or more photons. In the next section we introduce a 
direct probe of the Rydberg state population that will 
give more information about the details of this process. 

II. CASCADE FLUORESCENCE 

As a complementary diagnostic to trap loss observa- 
tions, we use a photon counting module and a narrow- 
band interference filter to detect the 420 nm decay pho- 
tons from the 6P 3 / 2 state to the 5Si/ 2 ground state. 
These photons are only observed coming from the MOT 
cloud under conditions of Rydberg excitation. The flu- 
orescent branching ratio to the 6P 3 / 2 from a high nS or 
nD Rydberg state is calculated to be about b r = 0.15, 
varying only slightly with principal quantum number. 
(The predominant channel is emission to the 5P 3 /2 state 
due to its much larger u> 3 factor in the emission rate, 
which easily compensates for a slightly smaller dipole 
matrix element. This channel is difficult to observe in 
the presence of strong MOT fluorescence at the same 
wavelength.) The probability of multiple photon cascade 
into the 6P3/2 state is small because the long- wavelength 
photons required are disfavored. P and F Rydberg states 
predominantly cascade into 5S, 6S, and 4D levels, all of 
which lie energetically below the 6P state. Even for states 
slightly above the 6P, such as the 5D, the predominant 
decay channel is to states below the 6P. Thus the cas- 
cade fluorescence is likely a reliable probe of the S or D 
Rydberg state populations, and is relatively insensitive 
to P or F Rydberg states. Figure shows the observed 
cascade counts as the excitation lasers are scanned across 
a Rydberg resonance. 

The observed cascade signal has a background, mainly 
from dark counts, of around 250/s and reaches a peak 



signal of about 10000/s for the 28 D excitation. We 
can compare these observed rates with our expected 
rates from the calculated excitation rates and fluores- 
cent branching ratios. Accounting for the finite collection 
solid angle £1 = 3x 10~ 3 , r\ = 3.4 % detection efficiency of 
the photomultiplier tube, and a calculated 6P 3 /2 - 5S 1 / 2 
emission branching ratio b§ = 0.31 we expect to observe 
a cascade count rate 

c 6 = R 2 Ngb r b 6 TiSl - Ar . = 18, 000/s (3) 

A r + Abb 

at the peak of the 28D excitation resonance. The ratio of 
spontaneous to total decay rates accounts for the effects 
of black-body radiation causing radiative transitions to 
nearby Rydberg states. These states are assumed, by 
the argument above, not to result in detected cascade 
photons. For the data seen in Fig. [2]we observe 10000/s, 
only 55% of the expected count rate. This suggests that 
roughly 1/2 of the Rydberg atoms are being transferred 
out of the excitation state by some other process. The 
ratio of expected to detected cascade counts for the other 
excitation states (30S, 43D, and 58D) are (0.5, 0.6, 0.6). 
Again, by the arguments from the previous section, the 
explanation for this cannot be inelastic Rydberg-Rydberg 
collisions, which would either produce an extremely large 
trap loss rate in contradiction to observations, or would 
necessarily vary greatly with principal quantum number. 

III. STIMULATED EMISSION PROBE 

For dipole blockade applications, one is particularly in- 
terested in the evolution of the blockaded Rydberg state. 
Furthermore, in the particular applications of dipole 
blockade to single atom and single photon sources |19j. 
stimulated emission is used to couple the blockaded Ry- 
dberg level to an intermediate atomic level. Thus it is 
natural for us to pursue the development of a stimu- 
lated emission probe of Rydberg dynamics. The inten- 
sity dependence of the signals produced by the stimulated 
emission probe allows further information to be obtained 
about the population dynamics. Such a probe is non- 
destructive, and has inherently high spectral resolution. 
While being less general and less sensitive than field ion- 
ization, it does have the potential to be applied to a num- 
ber of different states, subject to dipole selection rules. 

As shown in Fig. [TJ we apply a tunable diode laser in 
the range of 1013 to 1027 nm to perform stimulated emis- 
sion probing of the Rydberg states produced by the two- 
photon excitation. This "state probe" de-excites atoms 
from the Rydberg states to the 6P3/2 state, which sub- 
sequently decays via a 420 nm photon to the 5S ground 
state. The 6P 3 / 2 state is either detected in the same 
manner as the cascade light mentioned above, or by a 
reduction in trap loss (explained below) . As with the ex- 
citation lasers, the spatial size of the state-probe laser is 
comparable to the MOT cloud size. We deliver nearly 100 
mW of light to the atoms, corresponding to stimulated 
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FIG. 3: Stimulated emission probe scan across the 6P3/2 man- 
ifold. The cascade count background is around 1000 counts. 



FIG. 4: Loss rate dependence on stimulated emission probe 
intensity for the 28D state, showing short residence times for 
the Rydberg state produced by two-photon excitation. 



emission rates from 5xl0 5 to 5xl0 6 /s depending upon 
the Rydberg level involved. The natural decay rate of the 
6P3/2 level of 9.2xl0 6 /s is fast enough that population 
cannot accumulate there. The state-probe laser is stabi- 
lized to a 300 MHz optical spectrum analyzer that is itself 
locked to a 780 nm Rb saturated absorption resonance. 
It is tuned by changing the frequency of a double-passed 
acousto-optic modulator. 

When the stimulated emission probe is on resonance, 
atoms are returned to the ground state more quickly 
than they would otherwise spontaneously radiate from 
the long-lived Rydberg levels. This reduces the loss rate 
from the trap as the Rydberg atoms do not stay excited 
long enough for loss mechanisms such as black-body ion- 
ization and inelastic collisions to remove many of them 
from the trap. Keeping the excitation beams on reso- 
nance with the two-photon excitation, we can scan the 
frequency of the stimulated emission probe to observe 
the 6P3/2 hyperfrne manifold. Such a scan is shown in 
Fig.[3J When tuned on resonance with the 6P 3 / 2 F=3 hy- 
perfrne state, the MOT loss rate is reduced from 0.8/s to 
0.2/s, with the amount of reduction depending on probe 
intensity. 

The state-probe laser produces a reduction in the trap 
loss rates. An example of this for the 28D 5 / 2 Rydberg 
state is shown in Fig. 0J There the trap loss is shown as 
a function of the stimulated emission rate from the probe 
beam. The figure clearly demonstrates that by using a 
sufficiently high state-probe intensity the atoms can be 
returned to the ground state before ionization or inelastic 
collisions can occur, thus reducing the loss rate of the 
MOT. The surprising feature of the data is that much 
higher intensities are needed than would be expected if 
the dominant population transfer from the Rydberg level 
were spontaneous decay or black-body transfer (rates of 
4xl0 4 /s and 2xl0 4 /s, respectively) |32j| . The residence 
time of the Rydberg atoms is deduced to be nearly 10 
/is. 

Alternatively, we can measure the state probe laser ef- 
fects by detecting the number of 6P 3 / 2 decay photons 
from stimulated emission from the Rydberg state. The 
count rate data, scaled by the number of ground state 




0.0 0.5 1.0 1.5 2.0x10 

Stimulated Emission Rate (1/s) 



FIG. 5: Loss rate dependence on stimulated emission probe 
intensity for the 28D state, with the MOT magnetic field 
switched off during excitation. 

MOT atoms, is shown in Fig. M for the 28D state. The 
cascade and background signal have been subtracted for 
this data, thus there are no counts when the stimulated 
emission rate is zero. As with the repletion data of 
Fig. [U the counts saturate at a stimulated emission rate 
of around lxl0 5 /s, which implies that the other rates 
out of the excitation state must be on this order. 
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FIG. 6: Dependence of 6P3/2 decay counts on stimulated 
emission probe intensity for the 28D state. 

From the trap loss or state probe data, the residence 
time of the Rydberg state varies from about 10 /is at n = 
28 to about 50 /is at n — 58. Thus the population is being 
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transferred out of the initial Rydberg state faster than 
can be accounted for by spontaneous decay and black- 
body transfer (decay times of 25 /us and 50 /us at n = 28, 
respectively). The transfer times increase with n, just 
the opposite dependence as would be expected for energy 
transfer collisions between Rydberg states. These would 
be expected to decrease with n due to rapid increase in 
van der Waals interaction strengths. This counter trend 
thus implies that Rydberg-Rydberg collisions are not the 
process responsible for the fast transfer out of the initial 
Rydberg state. 

Interestingly, the 30S state has the highest transfer- 
rate to other Rydberg states. This again implies that 
atom-atom interactions do not dominate the process, as 
the strengths of van der Waals interactions tend to be 
larger for D states than for S states. 

We have checked that the transfer rates increase with 
increasing excitation rate. This confirms that some trans- 
fer process besides black-body radiation is occurring. 

As can be seen from Fig. [H there is a non-zero loss 
rate at high probe intensities where the stimulated emis- 
sion laser should fully deplete the original excitation 
state. There are several possible explanations for this. 
If there were a loss mechanism that was being enhanced 
by the state-probe laser, this would produce such a be- 
havior. However, processes such as photoionization or 
light-induced 6P-5S collisions can be estimated to be far 
too weak to account for this effect. Another possibility 
has to do with Zeeman precession in the nd states caus- 
ing population to accumulate in inaccessible magnetic 
sublevels. The magnetic field gradient used to confine 
the MOT atoms is large enough that precession between 
magnetic sublevels of the Rydberg atoms at the edges of 
the MOT occurs at a rate of several MHz, which is on 
the order of the stimulated emission rate. The result of 
this precession is that a fraction of the atoms - up to 1/3 
- move to a state that is inaccessible to the state-probe 
laser because of the dipole selection rules determined by 
the polarization of the state-probe laser. For linearly po- 
larized state-probe light, m = ±5/2 Zeeman levels can- 
not be excited to the 6P 3 / 2 state. Thus population that 
accumulates in these levels cannot be de-excited by the 
state-probe laser. 

This effect was verified by repeating the experiment 
with the MOT magnetic field being switched off for 10 
ms intervals and only switching the Rydberg excitation 
lasers during the times the field was off. This data, shown 
in Fig. [5] for n = 28 shows that Rydberg populations 
at high state probe intensities are markedly reduced as 
compared to when the magnetic field is on. The reduced 
signal-to-noise for this experiment made it possible to do 
this only for n = 28. 



IV. MODEL OF RYDBERG POPULATION 
DYNAMICS 

In the previous sections, we have described the basic 
processes that are evidently at work under the conditions 
of our experiment, and their experimental signatures. To 
further analyze the results, we present here a simplified 
model of the Rydberg dynamics and use it to extract the 
values of a few simple parameters from the data. Con- 
firming the interpretation in the previous section, we find 
that some process that does not cause trap loss neverthe- 
less transfers population out of the excitation Rydberg 
state on a time scale short compared to spontaneous or 
black-body lifetimes. The lengthening of this time scale 
with principal quantum number leads us to believe that 
it is not inelastic collisions between Rydberg atoms. 

The processes included in the model are illustrated in 
Fig. [71 We describe the system with a three state model: 
A ground state \g) with N g atoms, the excitation Ryd- 
berg state |r) with N r atoms, and an additional effective 
Rydberg state \s) that accounts for other states that are 
populated from state |r). 

The dynamics of the excitation state \r) depend 
on laser excitation and de-excitation, spontaneous and 
black-body radiation, and transfer to the other Rydberg 
states |s). Population enters |r) by excitation from the 
ground state at a rate R,2N g where i?2 is calculated as 
described in Section |TJ Spontaneous decay to low-lying 
levels occurs at a rate A r N r . Black-body radiation and 
other potential processes that transfer atoms to other 
Rydberg states occur at a rate jN r . There is also the 
possibility of trap loss (through ionization, for example) 
at a rate T r directly from state \r). Finally, de-excitation 
from the state-probe laser occurs at a rate R,3N r . Thus 

= R 2 N g - A r N r - R 3 N r - jN r - T r N r (4) 

is the rate equation for the excitation state population. 

The other Rydberg states are produced by collisional 
or radiative transfer from state \r) at the rate jN r and 
have an effective radiative lifetime A s . We also assume 
that these states can cause trap loss at a rate T s due 
to black-body ionization and other collisional processes. 
Thus they obey 

^ = 7 AT r - T S N S - A S N S (5) 
at 

We are assuming that transfer from s) back to \r) is 
unlikely. 

In addition to the radiative de-excitation and excita- 
tion processes with the Rydberg levels, the ground state 
population N g is affected by MOT loading (L) and loss 
(ro) processes that we assume are not materially changed 
when the Rydberg excitation lasers are on. The resulting 
rate equation for the ground state population is 

^JL = L- T N g - R 2 N g + (A r + R 3 )N r + A S N S (6) 
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FIG. 7: Diagram of the simplified model of Rydberg popula- 
tion dynamics. 



Plugging in the steady-state solutions to Eq.|3]and[S]gives 

(7) 



The loss-rate from the trap is increased by an amount 



A r + R 3 + 7 
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(8) 



This is a product of two factors. The first essentially 
measures the excitation rate of Rydberg atoms, modified 
by the de-excitation due to the state-probe laser. The 
second factor is the probability that the excited Rydberg 
atoms experience trap loss as opposed to radiatively de- 
caying back to the ground state. 

The model similarly predicts the count rate of 420 nm 
photons produced by the state-probe laser: 



(9) 



Thus the state-probe-induced 420 nm count rate can also 
be used to determine the transfer rate 7, with the data 
and fit for 28D shown in Fig. [6l 

The results of fitting our experimental data to this 
model are summarized in Table HI We list there the de- 
duced values of the primary unknowns 7 and T s , as well 
as the assumed values for the input parameters A r and 
A s . We note that both trap loss and 420 nm count rates 
can be used to extract 7, an important internal consis- 
tency check. The primary results are these: 

• The atoms transfer out of the excitation Rydberg 
state at a rate 7 that is substantially faster than 
spontaneous decay or black-body transfer rates. 

• The mechanism for population transfer, to an ex- 
cellent approximation, does not cause trap loss. 

• The population transfer rate decreases slowly with 
increasing principle quantum number, as opposed 
to the expected rapid increase if near-resonant en- 
ergy transfer collisions were the relevant mecha- 
nism. 



• The probability of trap loss is very small; most 
Rydberg excitations result in radiative repopula- 
tion of the ground state without trap loss. 

These conclusions from fitting the experimental data to 
the model are consistent with the simplified analyses pre- 
sented in Sections [11 11111 In the next two sections we will 
discuss their implications. 



TABLE I: A summary of transfer rate data using both count 
rate and loss rate methods. As a basis for comparison, the 
black-body transfer rate Abb, spontaneous emission rate A r , 
the estimated mean spontaneous emission rate A s from other 
Rydberg states, and the inferred trap loss rate F 3 are also 
included. All rates are in units of s _1 . 



State 7(counts) 


7(loss) Abb 




A s 


r s 


28D 1.2xl0 5 


1.3xl0 5 2.6xl0 4 


4.1 x 10 4 


3.1 x 10 4 


265 


43D 7.4 xlO 4 


7.2xl0 4 l.lxlO 4 


l.lxlO 4 


2.0 x 10 4 


602 


58D 2.6xl0 4 


2.0xl0 4 6.1xl0 3 


4.8xl0 3 


7.4 x 10 3 


433 


30S 3.9xl0 5 


5.0xl0 5 2.3xl0 4 


4.4xl0 4 


3.3 x 10 4 


83 



V. SUPERRADIANT POPULATION 
TRANSFER 

We have argued above that the transfer mechanism 
responsible for population transfer out of the excita- 
tion state in a few microseconds cannot be due to near- 
resonant energy transfer collisions between excited Ry- 
dberg atoms. L-changing collisions between Rydberg 
atoms and free electrons are another possible mechanism. 
Such collisions are unlikely to be energetic enough for the 
recoil to eject atoms from the trap and so would agree 
with much of the observed behavior. However, at our low 
excitation rates it is unlikely that electrons are present 
in sufficient quantities to cause the observed state trans- 
fer. With black-body ionization rates of 400/s (see Sec- 
tion [VlJ acting on a population of 10 4 Rydberg atoms, 
free electrons are produced at a rate of 4xl0 6 /s. Free 
electrons are typically fast-moving and dissipate at a rate 
of ~10 4 /s (33|, giving an average free electron population 
of ^400, which would not have a large impact on the 10 4 
atoms in the excitation state. We estimate that there 
is insufficient ionization to cause an electron trap [33j . 
In addition, the rate of l-changing collisions should in- 
crease with principle quantum number, in contrast with 
our observations. 

Having argued against collisional phenomena being re- 
sponsible for Rydberg energy transfer, we need a radia- 
tive mechanism to explain our results. Recently, Wang 
et al. [2c| observed superradiance in the measured life- 
times of Rydberg states. That superradiance could play 
an important role in Rydberg population dynamics can 
be understood by the following arguments. For principal 
quantum numbers > 20, the size of the MOT is less than 
the wavelength for radiative emission from the excita- 
tion state to nearby dipole-allowed states. If there are N 
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atoms initially in a particular Rydberg state, the collec- 
tive dipole moment is enhanced by a factor of N (assum- 
ing the atom cloud is much smaller than the wavelength 
of the emitted light). The emission rate is enhanced by 
a factor of N 2 , or by a factor of N on a per-atom basis. 
Since N is on the order of 10 4 for our experiment, the 
superradiant decay rate will be on the order of 



2 uj 3 Nd 2 

3 Tic 3 



4N 3 Ry 
3n 5 



a 



(10) 



which is 3 x 10 5 /s for 10 4 n — 50 Rydberg atoms. The 
low emission frequency is compensated for by the large 
number of cooperatively radiating atoms. Since the spon- 
taneous decay rate is proportional to n~ 3 , the relative 
importance of superradiance and spontaneous decay goes 
only as n -2 , a factor of only 4 for this experiment. 

Wang et al. [28[ developed a sophisticated theory of 
superradiance with application to cold Rydberg gases. 
Here we present a simplified model of superradiance to 
use as a interpretive guide. We base this model on Dicke's 
original work [34[ as elucidated by Gross and Haroche [H[ 
and Rehler and Eberly [36| . 

In considering the superradiant decay of an initial state 
|e) to a lower energy state \l), the Dicke approach intro- 
duces an effective collective spin state of the 2N-level 
system as \ JM), with iV e = J + M atoms in state |e) 
and Ni = J — M atoms in state \l). The radiation rate is 
found to be T el (J(J + 1) - M(M - 1)) = T e/ iV e (iV ; + 1). 
To extend the two- level case to our multi-level case, we 
assume that we can model superradiance with a set of 
rate equations 



dt 



J2 relN e (N t + 1) + Tl'MNe + 1) (11) 



Ke 



V>e 



This model reproduces the key features of superradiance: 
a large initial inversion radiates at r e ; (per atom) at 
first, then as the inversion is reduced the rate acceler- 
ates to a maximum of NY' el /A, occurring over a time 
(In N)/(T el N) < l/r e! 3d]. As an additional check on 
the model, we have simulated the experiment of Ref. [37| 
(which was performed at high temperatures and much 
smaller n) and our model reproduces the dominant fea- 
tures of the data shown there. 

The emission rates T e i are the rates for spontaneous 
emission multiplied by a cooperativity parameter C e i : 



2e 2 co 3 el (2Ji + 1) 
mc 3 (2J e + l) 



feu 



(12) 



where the f e i are the calculated absorption oscillator 
strengths. For a system of atoms whose spatial extent 
is on the order of the wavelength of the transition, su- 
perradiance will occur at a reduced rate. The wavelength 
of the transition to the nearest lower lying Rydberg state 
varies from 0.17 cm for 30S-29P to 2.8 cm for the 58D- 
59P transition. For low n, this is quite close to the 1 mm 
spatial extent of the MOT, and lower levels will have even 



shorter transition wavelengths. Following Ref. [36(, the 
cooperativity parameter for a uniform density system of 
N atoms in a volume V radiating in direction k is (in the 
N > 1 limit) 



1 

V 2 

9(sm(k e iR) - k e iRcos(k e iR)) 2 
(keiRf 



d 3 X / d 3 X j{k-ki)-{S-3) ( 13 ) 



(14) 



for a spherical uniform density cloud, where kR is the 
product of the wavenumber and the radius of the atomic 
sample. The cooperativity parameter is 1 for R <C A and 
decreases to for R » A. In practice, the cooperativity 
parameter becomes small around n = 20 for a 1 mm 
MOT. 

To account for superradiance, we replace the level \s) 
of the model of Section IIVI with a set of levels near in 
energy to \r) that are coupled to each other and to \r) by 
black-body radiation and by superradiance from Eq. 111! 
We find steady-state solutions of the resulting non-linear 
equations and from them deduce the effective Rydberg- 
Rydberg transfer rate from Eq. [TT]with e = r. 

Table II compares the observed Rydberg transfer rates 
and those predicted by our simulation. These rates do 
not drop off as quickly as would be expected from Equa- 
tion [lOl This is a result of the spatial factor approaching 
unity for the higher n-levels, balancing out the decrease in 
natural emission rate. Additionally as a result of this in- 
creasing spatial factor, our simulation indicates that the 
atoms are transferred predominantly to f-states for the 
43d and 58d levels. This transfer to f-states could explain 
the higher total loss rate from the trap for these levels. 
We emphasize that there are no adjustable parameters in 
our simplified superradiance model, and yet it naturally 
predicts the order of magnitude of the Rydberg- Rydberg 
transfer rates. A more sophisticated model would be ex- 
pected to explain the variation seen. 



TABLE II: A comparison of Rydberg-Rydberg transfer rates 
deduced by comparing measurements deduced from the model 
of Section IIVI with those predicted by the superradiance 
model. All units are s" 1 . 



State 7(calculated) 7(expt) 



28D 


1.7xl0 5 


1.3xl0 5 


43D 


2.4x10 s 


7.4xl0 4 


58D 


1.2xl0 5 


2.0xl0 4 


30S 


2.2xl0 5 


5.0xl0 5 



VI. BLACK-BODY IONIZATION 

So far we have not discussed the actual mechanism for 
trap loss. The deduced trap loss rates are quite modest, 
typically 400/s on a per Rydberg atom basis. This is 
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close to what would be expected from black-body ion- 
ization. Black-body ionization rates were recently cal- 
culated in Ref. 38] for the various Rydberg levels. As- 
suming Rydberg n-level distributions from our superradi- 
ance model, we find expected photoionization rates from 
black-body radiation to be as shown in Table ITIT1 For the 
D states, there is reasonable agreement between our de- 
duced experimental rates and the predicted values, sug- 
gesting that black-body ionization probably composes a 
large portion of the total loss from the trap. For the s- 
state, the predicted trap loss is greater than observed by 
about a factor of 3, for which we have no explanation. 

TABLE III: Comparison of inferred Rydberg trap loss rates 
with black-body ionization rates. Tbbi is the rate of black- 
body ionization for each of these states calculated from 



are in 


units of 


s" 1 . 


State T s calc 




28D 


212 


322 


43D 


470 


720 


58D 


329 


457 


30S 


77 


265 



VII. DISCUSSION 

In this paper, we have presented results on the dy- 
namics of low-density ultracold Rydberg clouds excited 
using two-photon absorption from a MOT. By looking 
at trap loss, radiative cascade, and state-selective stim- 
ulated emission we find that the probability of the Ry- 
dberg atoms undergoing collisional loss before decaying 
back to the ground state is small. Additionally, the Ryd- 
berg state produced by the two-photon excitation is de- 
pleted by some mechanism on a time scale significantly 
shorter than can be explained by black-body transfer or 
single-atom radiative decay. This time scale increases 
with increasing principal quantum number, in contradic- 



tion to expectations if inelastic Rydberg-Rydberg colli- 
sions were responsible for the rapid state transfer. On 
the other hand, estimates of collective superradiant light 
emission and a simplified model thereof suggest that this 
mechanism can explain the rapid population transfer ob- 
served in this experiment. The overall trap loss rates are 
consistent with expectations from black-body ionization. 

Before the experiment of Wang et a/.[28| the effects 
of superradiance were not appreciated for the popula- 
tion dynamics of high density ultracold Rydberg atoms. 
Superradiance can happen on very rapid time scales, es- 
pecially under strong Rydberg excitation as achieved in 
a number of previous experiments. In fact, superradi- 
ance can easily occur on sub-microsecond time scales. 
A key consequence of superradiance is that it populates 
Rydberg states lying energetically below the state be- 
ing excited by the laser. Indeed, in the experiments of 
Ref. [l3[ and Ref. [3f| , population clearly moves to lower 
lying Rydberg levels on a fast time scale, consistent with 
the hypothesis of superradiant transfer. 

An additional consequence of fast superradiant popu- 
lation transfer is that it provides a mechanism for rapid 
population of states of neighboring orbital angular mo- 
mentum I from the excitation state. Pairs of atoms 
with 61 = ±1 interact at long range via the i? -3 res- 
onant dipole-dipole interaction, not the usual R~ 6 van 
der Waals interaction. The much stronger collision inter- 
actions between these atoms may explain the very rapid 
time scales for plasma formation in a number of experi- 
ments where resonant dipole-dipole interactions were not 
purposely produced using external fields. 



Acknowledgments 

J. Strabley contributed to the early stages of this work. 
We appreciate helpful discussions with M. Saffman and 
D. Yavuz. This work was supported by the National 
Science Foundation and NASA. 



[1] I. Mourachko, D. Comparat, F. de Tomasi, A. Fioretti, 
P. Nosbaum, V. M. Akulin, and P. Pillet, Phys. Rev. 
Lett. 80, 253 (1998). 

[2] W. R. Anderson, J. R. Veale, and T. F. Gallagher, Phys. 
Rev. Lett. 80, 249 (1998). 

[3] S. Farooqi, D. Tong, S. Krishnan, J. Stanojevic, 
Y. Zhang, J. Ensher, A. Estrin, C. Boisseau, R. Cote, 
E. Eyler, et al., Phys. Rev. Lett. 91, 183002 (2003). 

[4] K. Afrousheh, P. Bohlouli-Zanjani, D. Vagale, A. Mug- 
ford, M. Fedorov, and J. D. D. Martin, Phys. Rev. Lett. 
93, 233001 (2004). 

[5] K. Afrousheh, P. Bohlouli-Zanjani, J. D. Carter, A. Mug- 
ford, and J. D. D. Martin, Phys. Rev. A 73, 063403 
(2006). 

[6] T. Vogt, M. Viteau, J. Zhao, A. Chotia, D. Comparat, 



and P. Pillet, Phys. Rev. Lett. 97, 083003 (2006). 
[7] P. Bohlouli-Zanjani, J. A. Petrus, and J. D. D. Martin, 

Phys. Rev. Lett. 98, 203005 (2007). 
[8] W. Li, P. J. Tanner, and T. F. Gallagher, Phys. Rev. 

Lett. 94, 173001 (2005). 
[9] K. Singer, M. Reetz-Lamour, T. Amthor, S. Foiling, 

M. Tscherneck, and M. Weidemiiller, J. Phys. B 38, S321 

(2005). 

[10] T. Amthor, M. Reetz-Lamour, S. Westermann, J. Den- 
skat, and M. Weidemiiller, Phys. Rev. Lett. 98, 023004 
(2007). 

[11] T. J. Carroll, S. Sunder, and M. W. Noel, Phys. Rev. A 
73, 032725 (2006). 

[12] M. P. Robinson, B. L. Tolra, M. W. Noel, T. F. Gal- 
lagher, and P. Pillet, Phys. Rev. Lett. 85, 4466 (2000). 



9 



[13] A. Walz-Flannigan, J. R. Guest, J. H. Choi, and 

G. Raithel, Phys. Rev. A 69, 063405 (2004). 
[14] D. Jaksch, J. I. Cirac, P. Zoller, S. L. Rolston, R. Cote, 

and M. D. Lukin, Phys. Rev. Lett. 85, 2208 (2000). 
[15] M. D. Lukin, M. Fleischhauer, R. Cote, L. M. Duan, 

D. Jaksch, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 87, 

037901 (2001). 

[16] J. Sebby-Strabley, R. T. R. Newell, J. O. Day, E. Brekke, 

and T. G. Walker, Phys. Rev. A 71, 021401(R) (2005). 
[17] D. D. Yavuz, P. B. Kulatunga, E. Urban, T. A. Johnson, 

N. Proite, T. Henage, T. G. Walker, and M. Saffman, 

Phys. Rev. Lett. 96, 063001 (2006). 
[18] E. Brion, K. Molmer, and M. Saffman, Phys. Rev. Lett. 

99, 260501 (pages 4) (2007). 
[19] M. Saffman and T. G. Walker, Phys. Rev. A 66, 065403 

(2002). 

[20] M. Saffman and T. G. Walker, Phys. Rev. A 72, 042302 
(2005). 

[21] T. G. Walker and M. Saffman, arXiv:0712.3438vl [quant- 
ph] (2008). 

[22] R. Heidemann, U. Raitzsch, V. Bendkowsky, B. Butscher, 
R. Low, L. Santos, and T. Pfau, Phys. Rev. Lett. 99, 
163601 (2007). 

[23] K. Singer, M. Reetz-Lamour, T. Amthor, L. G. Marcassa, 

and M. Weidemiiller, Phys. Rev. Lett. 93, 163001 (2004). 
[24] D. Tong, S. M. Farooqi, J. Stanojevic, S. Krishnan, Y. P. 

Zhang, R. Cote, E. E. Eyler, and P. L. Gould, Phys. Rev. 

Lett. 93, 063001 (2004). 
[25] T. C. Liebisch, A. Reinhard, P. R. Berman, and 

G. Raithel, Phys. Rev. Lett. 95, 253002 (2005), Phys. 

Rev. Lett. 98, 109903(E) (2007). 
[26] T. A. Johnson, E. Urban, T. Henage, L. Isenhower, D. D. 

Yavuz, T. G. Walker, and M. Saffman, arXiv:0711.0401 



(2008). 

[27] R. Heidemann, U. Raitzsch, V. Bendkowsky, B. Butscher, 
R. Low, and T. Pfau, Phys. Rev. Lett. 100, 033601 
(2008). 

[28] T. Wang, S. F. Yelin, R. Cote, E. E. Eyler, S. M. Farooqi, 
P. L. Gould, M. Kostrun, D. Tong, and D. Vrinceanu, 
Phys. Rev. A 75, 033802 (pages 4) (2007). 

[29] S. K. Dutta, D. Feldbaum, A. Walz-Flannigan, J. R. 
Guest, and G. Raithel, Phys. Rev. Lett. 86, 3993 (2001). 

[30] A. Reinhard, T. C. Liebisch, B. Knuffman, and 
G. Raithel, Phys. Rev. A 75, 032712 (2007). 

[31] B. K. Teo, D. Feldbaum, T. Cubel, J. R. Guest, P. R. 
Berman, and G. Raithel, Phys. Rev. A 68, 053407 (2003). 

[32] T. F. Gallagher, Rydberg Atoms (Cambridge University 
Press, 1994). 

[33] T. C. Killian, S. Kulin, S. D. Bergeson, L. A. Orozco, 
C. Orzel, and S. L. Rolston, Phys. Rev. Lett. 83, 4776 
(1999). 

[34] R. H. Dicke, Phys. Rev. 93, 99 (1954). 
[35] M. Gross and S. Haroche, Phys. Rep. 93 (1982). 
[36] N. E. Rehler and J. H. Eberly, Phys. Rev. A 3, 1735 
(1971). 

[37] P. R. F. F Gounand, M Hugon and 
J. Berlande, Journal of Physics B: Atomic 
and Molecular Physics 12, 547 (19 79), URL 
|http : //stacks . iop . org/0022-3700/12/547| 

[38] I. I. Beterov, D. B. Tretyakov, I. I. Ryabtsev, A. Ekers, 
and N. N. Bezuglov, arXiv.org:physics/0702192 (2007). 

[39] W. Li, M. W. Noel, M. P. Robinson, P. J. Tanner, T. F. 
Gallagher, D. Comparat, B. Laburthe Tolra, N. Van- 
haecke, T. Vogt, N. Zahzam, et al., Phys. Rev. A 70, 
042713 (2004). 



